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EXECUTIVE  SUMMARY 


The  prediction  of  electromagnetic  field  strengths  at  large  distances  from 
a  radiating  antenna  is  important  to  designers  and  users  of  communications 
systems.  Such  field  strengths  can  be  greatly  influenced  by  the  variation  of 
the  modified  atmospheric  refractivity  M  with  height  z.  If  the  gradient  of  M 
is  negative  over  a  range  of  z,  an  extreme  enhancement  in  field  strength  may  be 
realized  at  large  distances  from  the  source.  The  layer  of  the  atmosphere 
responsible  for  the  enhancement  is  known  as  a  duct.  When  the  ground  is  the 
lower  boundary  of  the  ducting  layer,  the  duct  is  called  a  surface,  or  ground- 
based,  duct;  otherwise,  it  is  denoted  as  an  elevated  duct.  One  common  type  of 
surface  duct  is  found  over  the  ocean  and  is  known  as  an  evaporation  duct. 

Since  the  propagation  of  electromagnetic  waves  is  sensitive  to  the  phase  of 
plane  wave  spectral  components,  which  can  be  greatly  influenced  by  the  form  of 
the  reflecting  ground,  ground  roughness  is  expected  to  have  an  important 
effect  on  the  field  strengths. 

As  a  first  step  toward  developing  a  mathematical  model  for  duct 
propagation  over  a  general  rough  earth,  the  problem  considered  herein  is  of  a 
perfectly  conducting  periodic  ground  with  a  small  roughness  amplitude  on  the 
lower  boundary  of  a  duct.  Specifically,  the  wave  equation  is  developed  and 
solved  for  the  two-dimensional  case  of  a  long  horizontal  source.  This  method 
is  applied  to  the  problem  of  a  long  horizontal  source  radiating  over  a 
sinusoidally  shaped,  perfectly  conducting  ground  within  a  surface  duct 
environment.  The  eigenvalues  of  the  "rough-walled  waveguide"  are  found  by 
expanding  a  characteristic  determinant  about  zero  roughness  amplitude  to 
obtain  an  "equivalent  zero  mode  matrix".  It  is  shown  that  for  the  specular 
component  of  the  reflection,  the  problem  can  be  considered  in  terms  of  an 
effective  reflection  coefficient  that  accounts  for  re-scattering  of  energy 
back  into  the  specular  direction. 
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An  eigenvalue  perturbation  method  is  developed  using  the  equivalent 
matrix  that  permits  the  calculation  of  eigenvalues  for  the  rough-earth  case 
using  those  of  the  smooth-earth  case.  It  is  shown  that  the  eigenvalue 
perturbation  scheme  is  as  accurate  as  directly  obtaining  the  eigenvalues  from 
the  equivalent  matrix  but  requires  much  less  computation  time.  The 
relationship  between  the  two  methods  has  been  detailed,  and  the  eigenvalue 
perturbations,  as  well  as  roughness  loss  rates,  have  been  shown  to  be 
consistent  with  results  obtained  using  a  rough-earth  form  of  the  fundamental 
waveguide  mode  equation. 

The  methods  documented  in  this  report  are  a  first  step  toward  developing 
a  mathematical  model  of  propagation  under  the  given  conditions.  They  can  be 
considered  a  research  model  applicable  to  a  limited  set  of  circumstances . 
Further  work  is  needed  to  extend  the  model  and  make  it  more  widely  applicable. 
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PREFACE 


The  Electromagnetic  Compatibility  Analysis  Center  (ECAC),  a  Department  of 
Defense  facility,  was  established  to  provide  advice  and  assistance  on 
electromagnetic  compatibility  matters  to  the  Secretary  of  Defense,  the  Joint 
Chiefs  of  Staff,  the  military  departments,  and  other  DoD  components.  The 
Center,  located  at  North  Severn,  Annapolis,  Maryland  21402,  is  under  the 
policy  control  of  the  Assistant  Secretary  of  Defense  for  Communication, 
Command,  Control,  and  Intelligence,  and  the  Chairman,  Joint  Chiefs  of  Staff, 
or  their  designees,  who  jointly  provide  policy  guidance,  assign  projects,  and 
establish  priorities.  ECAC  functions  under  the  executive  direction  of  the 
Secretary  of  the  Air  Force,  and  the  management  and  technical  direction  of  the 
Center  are  provided  by  military  and  civil  service  personnel.  The  technical 
support  function  is  provided  through  an  Air  Force-sponsored  contract  with  the 
I IT  Research  Institute  (IITRI). 

To  the  extent  possible,  all  abbreviations  and  symbols  used  in  this  report 
are  taken  from  the  American  National  Standards  Instutute,  Inc.,  American 
National  Standard  ANSI  (Y10.19)  1969  letter  Symbols  for  Units  Used  in  Science 
and  Technology. 

Users  of  this  report  are  invited  to  submit  comments  that  would  be  useful 
in  revising  or  adding  to  this  material  to  the  Director,  ECAC,  North  Severn, 
Annapolis,  Maryland  21402-1187,  Attention:  XM. 
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duct 


elevated  duct 


evaporation  duct 


Floquet  series 


index  of  refraction,  n 


modified  refractivity ,  M 


perturbation  expansion 


GLOSSARY 

-  a  nonstandard  condition  in  the  troposphere  whereby 
a  layer  is  formed,  the  refractivity  gradient  of 
which  is  less  than  -157  N  units/km,  causing  the 
ray  curvature  of  a  radio  wave  to  be  greater  than 
the  earth's  curvature,  with  the  result  that  the 
wave  is  trapped  in  the  layer  and  propagated  beyond 
its  normal  range 

-  a  duct,  the  bottom  of  which  is  above  the  earth's 
surface 

-  a  surface  duct  formed  above  a  body  of  water  when 
evaporation  causes  a  moist  air  layer  to  form  under 
a  dry  layer 

-  a  discrete  sum  of  plane  waves  used  as  part  of  the 
solution  of  a  differential  equation  with  periodic 
boundary  conditions 

-  ratio  of  the  phase  velocity  of  an  electromagnetic 
wave  in  free  space  to  the  phase  velocity  in  a 
medium 

-  defined  by  the  equation  M  =  N  +  (h/r)'10  ,  where  N 
is  the  refractivity,  h  is  the  height,  and  r  is  the 
earth's  radius.  M  is  defined  in  terms  of  N  so 
that  its  gradient  is  negative  in  a  duct. 

-  expansion  of  the  rough-earth  eigenvalues  in  terms 
of  the  smooth-earth  ones  as  a  series  in  the 
roughness  height  parameter,  valid  for  small  values 
of  the  roughness  height 

xi 
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GLOSSARY  (Continued) 

refractivity ,  N  -  defined  for  convenience  by  the  equation 

(n-l)'IO^,  where  n  is  the  index  of  refraction 

refractivity  profile  -  a  plot  of  refractivity  (either  N  or  M)  versus 

height,  where  refractivity  is  the  abscissa  and 
height  is  the  ordinate 

specular  reflection  -  reflection  that  is  the  same  type  as  that  caused  by 

smooth  surfaces  and  has  the  following 
properties:  it  is  directional  (angle  of  incidence 
equals  the  angle  of  reflection),  its  phase  is 
coherent,  and  its  fluctuations  have  a  relatively 
small  amplitude 

surface  duct  -  a  duct  with  a  bottom  that  is  on  the  earth's 

surface  and  a  modified  refractivity  gradient  that 
is  negative  from  the  earth's  surface  to  the  duct 
top 
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SECTION  1 
INTRODUCTION 


BACKGROUND 


The  prediction  of  electromagnetic  field  strengths  at  large  distances  from 

a  radiating  antenna  is  important  to  designers  and  users  of  communications 

systems.  Such  field  strengths  can  be  greatly  influenced  by  the  variation  of 

the  modified  atmospheric  refractivity  M  with  height  z.  If  dM/dz  is  negative 

over  a  range  of  z,  the  field  strength  can  be  enhanced  at  large  distances  from 

the  source.  The  layer  of  the  atmosphere  responsible  for  the  enhancement  is 

called  a  duct.1-1  When  the  ground  is  the  lower  boundary  of  the  ducting  layer, 

the  duct  is  called  a  surface,  or  ground-based,  duct.  Otherwise,  it  is  denoted 

as  an  elevated  duct.  One  common  type  of  surface  duct  is  found  over  the  ocean 

1  P 

and  is  known  as  an  evaporation  duct. 

Numerical  methods  for  predicting  field  strengths  at  large  distances  from  a 
radiating  source  in  surface  and  elevated  ducts  have  been  documented.1-^’1-11 
These  methods  are  based  on  the  waveguide  mode  theory  of  wave  propagation1--* 
and  assume  horizontal  homogeneity  of  the  tropospheric  layers  and  a  smooth 
earth.  Since  the  ducted  modes  are  sensitive  functions  of  the  phase  of  plane 


1  ^err,  D.  E.,  "Propagation  of  Short  Radio  Waves,"  MIT  Radiation  Laboratory 
Series,  Vol.  13,  McGraw  Hill,  New  York,  1951. 

1  p 

Rotheram,  S.,  "Radiowave  Propagation  in  the  Evaporation  Duct,"  The  Marconi 
Review,  Vol.  37,  No.  192,  1974,  pp.  18-40. 

1-^Marcus,  S.  W.,  "A  Model  to  Calculate  EM  Fields  in  Tropospheric  Duct 

Environments  at  Frequencies  through  SHF,"  Radio  Science,  Vol.  17,  No.  5, 

1982,  pp.  885-901. 

1-i1Marcus,  S.  W.  ,  and  Stuart,  W.  D.,  A  Model  to  Calculate  EM  Fields  in 

Tropospheric  Duct  Environments  at  Frequencies  through  SHF,  ESD  Technical 
Report  81-102,  AD-A107710,  DoD  ECAC,  Annapolis,  MD,  September  1981. 

1-^Budden,  K.  G.,  The  Wave  Guide  Mode  Theory  of  Wave  Propagation,  Logos 
Press,  Ltd.,  London,  1961. 
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wave  spectral  components,  which  can  be  greatly  influenced  by  the  form  of  the 
reflecting  ground,  ground  roughness  is  expected  to  have  an  important  effect  on 
these  modes. 

A  method  for  considering  ground  roughness  in  waveguide  models  of  duct 
propagation  has  been  developed  and  utilized  by  Rotheram  (Reference  1-2)  and 
Hitney.1-®  In  the  expression  for  the  boundary  conditions  at  the  ground,  this 
method  replaces  the  Fresnel  reflection  coefficient  by  the  effective  reflection 
coefficient  of  a  randomly  rough  surface  as  developed  by  Rice1-^  and  Ament1-® 
and  verified  by  Beard. 1-^ 

For  the  case  of  a  horizontally  polarized  wave  incident  on  a  perfect 
conductor,  the  Fresnel  reflection  coefficient  is  -1,  while  the  effective 
reflection  coefficient  is 

R  =  -1  +  1/2  (A<fr)2  (1-1) 

where  the  Rayleigh  factor  A$  is  defined  as 

A4>  =  2  k  h  sin  9  (1-2) 

k  is  the  free-space  wave  number,  h  is  the  root-mean-square  (rms)  bump  height 
of  the  roughness,  and  6  is  the  grazing  angle  of  the  ray  representing  the 
incident  plane  wave. 


1-®Hitney,  H. ,  et  al . ,  "Tropospheric  Radio  Propagation  Assessment," 

Proceedings  of  the  IEEE,  Vol.  73,  No.  2,  1985,  pp.  265-283. 

1-^Rice,  S.  0.,  "Reflection  of  Electromagnetic  Waves  from  Slightly  Rough 

Surfaces,"  The  Theory  of  Electromagnetic  Waves,  M.  Kline,  ed.,  Interscience 
Publishers,  New  York,  NY,  1951,  p.  351. 

1-®Ament,  W.  D.,  "Toward  a  Theory  of  Reflection  by  a  Rough  Surface," 
Proceedings  of  the  IRE,  Vol.  41,  No.  1,  1953,  pp.  142-146. 

1-^Beard,  C.  I.,  "Coherent  and  Incoherent  Scattering  of  Microwaves  from  the 
Ocean,"  IRE  Transactions  on  Antennas  and  Propagation,  Vol.  AP-9,  No.  5, 
1961,  pp.  470-483. 
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This  effective  coefficient  considers  only  reflection  in  the  specular 
direction.  It  is  an  expression  of  the  amount  of  energy  scattered  out  of  that 
direction  by  the  ground  roughness  and  it  contains  the  assumption  of  a  standard 
atmosphere  with  no  ducts.  This  method  of  considering  ground  roughness  when 
ducts  are  present  is  popular  in  spite  of  its  "semi-quantitativeness" 

(Reference  1-6).  One  difficulty  in  justifying  its  use  is  based  on  the 
dependency  of  waveguide  modes  upon  continuous  reflection  from  the  waveguide 
wails  (Reference  1-5).  Thus,  a  portion  of  the  energy  scattered  out  of  the 
specular  direction  at  one  reflection  from  the  ground  may  be  scattered  back 
into  the  original  specular  direction  at  a  later  reflection  from  the  ground. 

In  addition,  because  of  the  nature  of  the  duct,  energy  scattered  out  of  the 
specular  direction  may  continue  to  propagate  and  may  represent  a  mode  in  its 
own  right.  Finally,  the  effective  reflection  coefficient  was  derived  and 
experimentally  verified  for  real  angles  of  incidence.  Its  carryover  to  the 
complex  eigenangles  characteristic  of  the  waveguide  formulation  of  duct 
propagation  has  not  been  validated. 

OBJECTIVE 

The  objective  of  this  analysis  was  to  develop  a  first-degree 
approximation  for  estimating  the  electromagnetic  field  strengths  at  large 
distances  from  a  radiating  antenna  that  is  within  a  surface  duct  situated  over 
a  periodically  rough  ground. 

APPROACH 


As  a  first  step  toward  developing  a  mathematical  model  for  duct 
propagation  over  a  general  rough  earth,  the  problem  considered  herein  is  of  a 
perfectly  conducting  periodic  ground  with  a  small  roughness  amplitude  below  a 
duct  environment.  The  wave  equation  is  developed  and  solved  for  the  two- 
dimensional  case  of  a  long,  horizontal  line  source.  The  eigenvalues  of  the 
solution  are  the  roots  of  an  infinite  determinant  that  reduce  to  the  smooth- 
earth  eigenvalues  as  the  ground  roughness  approaches  zero.  Because  it  is 
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infinite-dimensional,  the  determinant  cannot  be  evaluated  directly  and  the 
eigenvalue  must  be  calculated  using  other  means. 

Towards  this,  an  expression  for  the  ground  reflection  coefficient 
analogous  to  Equation  1-1  is  derived  that  accounts  for  the  dispersive  nature 
of  the  ground-reflected  field,  as  well  as  the  multiple  reflections 
characteristic  of  waveguide  modes.  This  effective  reflection  coefficient  is 
obtained  for  the  case  of  horizontally  polarized  radio  waves  propagating  over  a 
sinusoidally  shaped,  perfectly  conducting  ground  within  a  surface  duct 
environment . 

This  reflection  coefficient  is  then  used  in  conjunction  with  the 
fundamental  waveguide  mode  equation  to  derive  an  eigenvalue  perturbation 
method  wherein  the  eigenvalues  for  propagation  over  rough  earth  can  be 
obtained  easily  from  the  smooth-earth  eigenvalues.  A  simple  numerical 
expression  is  obtained  for  this  eigenvalue  perturbation  in  terms  of  ground  and 
atmospheric  reflection  coeff icients .  Although  the  fundamental  waveguide  mode 
equation  approach  is  entirely  equivalent  to  the  formulation  based  on  the  wave 
equation,1-^  the  waveguide  mode  formulation  provides  greater  physical  insight 
into  the  waveguide  phenomena.  This  waveguide  mode  formulation  also  leads  to  a 
closed-form  expression  for  the  rate  of  modal  attenuation  due  to  roughness  for 
trapped  modes  in  a  surface  duct  formed  by  a  bilinear  refractivity  profile. 

Another  expression  for  the  eigenvalues  is  then  derived  by  expanding  a 
characteristic  determinant  about  zero  roughness  amplitude  to  obtain  an 
"equivalent  matrix"  for  propagation  over  rough  earth  which  is  identical  in 
form  to  the  characteristic  matrix  for  smooth-earth  propagation.  It  is  shown 
that  for  the  specular  component  of  the  reflection,  the  problem  can  be 
considered  in  terms  of  an  effective  reflection  coefficient  that  accounts  for 
rescattering  of  energy  back  into  the  specular  direction.  (The  determination 


^"^Marcus,  S.  W.,  "Propagation  in  a  Surface  Duct  Over  a  Two-Dimensional 
Sinusoidal  Earth,"  Radio  Science,  Vol.  23,  No.  6,  1988,  pp.  1039-1047. 


1-4 


ECAC-TR-88-001 


Section  1 


of  field  contributions  due  to  the  portion  of  energy  scattered  out  of  the 
specular  direction  and  not  rescattered  back  into  that  direction  will  be 
relegated  to  a  future  study.) 

A  perturbation  technique  is  used  to  derive  an  expression  for  the 
eigenvalues  of  the  rough-earth  case  in  terms  of  the  smooth-earth  eigenvalues 
from  the  equivalent  matrix.  This  eliminates  the  necessity  of  a  separate, 
time-consuming  search  for  the  roots  of  the  characteristic  determinant  of  the 
equivalent  matrix  for  each  set  of  roughness  parameters.  The  accuracy  of  this 
perturbation  method  is  discussed,  particularly  as  it  relates  to  the 
attenuation  of  trapped  modes  in  surface  ducts  due  to  the  roughness. 

Numerical  field  strength  predictions  are  presented  with  and  without  the 
sinusoidal  ground  roughness,  which  shows  that  attenuation  due  to  roughness  can 
be  significant.  A  comparison  is  also  made  between  predictions  for  which  the 
complete  reflection  coefficient  is  used  and  those  for  which  only  scattering 
out  of  the  specular  direction  is  considered. 

ORGANIZATION  OF  THE  DOCUMENT 


In  Section  2,  the  field  equation  for  a  horizontally  polarized  wave 
propagating  in  a  duct  situated  over  a  perfectly  conducting,  periodically  rough 
ground  with  small  amplitude  is  derived.  The  solution  to  this  equation  is 
given  in  the  form  of  a  Floquet  series.  The  formal  solution  is  dependent  on 
determining  eigenvalues  that  are  the  roots  of  an  infinite  determinant.  It 
remains  to  develop  practical  means  of  determining  the  roots.  Theoretical 
concepts  required  to  achieve  this  are  derived  in  Section  3  and  the  actual 
solution  is  obtained  in  Section  4. 

In  Section  3,  a  waveguide  mode  formulation  is  used  to  derive  an  effective 
ground  reflection  coefficient  when  the  roughness  is  sinusoidal.  This 
coefficient  is  used  in  the  fundamental  waveguide  mode  equation  to  determine 
eigenvalues  of  the  rough  waveguide  as  a  perturbation  of  the  smooth-earth 
values.  The  roughness  is  found  to  cause  attenuation  of  trapped  modes,  and 
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this  attenuation  is  expressed  in  a  simple  manner.  This  expression  reduces  to 
a  closed  form  when  the  duct  is  approximated  by  a  bilinear  refractivity 
profile. 

For  the  same  sinusoidal  roughness,  the  specular  contribution  to  the  field 
equations  of  Section  2  are  derived  in  Section  4.  This  is  done  by  deriving  an 
"equivalent  zero  mode  matrix"  from  the  infinite  determinant  using  the 
effective  reflection  coefficient.  This  accounts  for  scattering  both  out  of 
and  into  the  specular  direction.  A  computationally  efficient  means  of 
obtaining  the  relevant  eigenvalues  from  the  equivalent  matrix  as  perturbations 
of  the  smooth-earth  ones  are  derived.  It  should  be  noted  that  this 
perturbation  expansion  is  not  the  one  mentioned  in  the  previous  paragraph. 

Once  again,  expressions  for  loss  rates  are  obtained. 

In  Section  5,  numerical  results  obtained  directly  from  the  equivalent 
matrix  and  from  the  eigenvalue  perturbation  methods  are  compared.  This  is 
done  by  first  plotting  roughness  loss  rates  of  individual  trapped  modes, 
obtained  from  the  two  methods,  and  then  plotting  field  strengths  using  the 
four  lowest  order  modes. 

Section  6  contains  a  summary  of  the  report. 
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SECTION  2 

GENERAL  FORMULATION 


GENERAL 


References  1-3  and  1-4  contain  the  mathematical  formulation  for 
describing  an  electromagnetic  wave  propagation  in  a  tropospheric  duct  over  a 
smooth  earth.  This  formulation  will  now  be  extended  to  include  a  periodically 
rough  boundary  condition  between  the  earth  and  the  troposphere.  Solutions  to 
the  field  equations  will  be  obtained  using  a  Floquet  series.  These  solutions 
will  be  dependent  on  determining  the  eigenvalues  that  are  the  roots  of  an 
infinite  determinant.  Determination  of  the  eigenvalues  will  be  discussed  in 
Section  4. 

PROPAGATION  MEDIA  AND  ROUGH  BOUNDARY 


The  modified  index  of  refraction  m(z)  is  defined  as 

m(z)  =  M(z)  x  10"6  +  i  =  n(z)  +  z/a  (2-1) 

where  n(z)  is  the  actual  index  of  refraction;  a  is  the  radius  of  the  earth  in 
the  same  units  as  the  height  z;  and  the  z/a  term  compensates  for  earth 
curvature,  thereby  enabling  the  use  of  Cartesian  coordinates  (Reference  1-1). 

The  modified  refractivity  profile  (i.e.,  its  variation  with  height)  will 
be  approximated  by  a  continuous  piecewise  linear  profile  with  L  sections  (see 
Figure  2-1).  Each  section  will  represent  an  atmospheric  layer,  or  region,  the 
boundaries  of  which  are  parallel  to  a  flat  earth  located  at  z  =  Zq  =  0.  The 
interface  between  the  i-th  and  i-th  +  1  layer  is  located  at  z  : 

(1  <  i  <  L  -  1),  with  the  layer  i  =  1  closest  to  the  ground.  The  L-th  layer 
is  unbounded  in  the  positive  z  direction. 
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In  each  atmospheric  layer, 


2 .  . 
mi  (z) 


( z  -  H . )  tan  a. ,  1  <  i  <  L 

i  l  ’ 


(2-2) 


where  is  the  value  of  z  at  which  m^(z)  would  equal  unity  and  the  slope 
(tan  a^)/2  of  the  m ^  versus  z  curve  is  assumed  small. 

It  is  assumed  that  the  boundary  between  the  atmosphere  ana  the  infinitely 
conductive  ground  is  rough  in  a  periodic  manner  with  spatial  period  d.  This 
boundary  can  be  described  in  terms  of  a  Fourier  series. 

FIELD  EQUATIONS 


The  source  of  a  horizontally  polarized  wave  can  be  considered  to  be  a 
linear  density  of  magnetic  dipoles  j5  of  infinite  length  that  is  uniformly 
distributed  along  a  line  parallel  to  the  y  axis,  oriented  in  the  z  direction, 
and  located  at  x  =  x' ,  z  =  zT  (see  Figure  2-1).  In  a  laterally  homogeneous 
medium,  the  fields  due  to  such  a  source  may  be  obtained  from  the  z-directed 
magnetic  Hertz  potential  vector: 

n  =  n(x,z)  z  (2-3) 


and 


£  =  -ju)u0vxn,  fi  =  vxvxn  (2-4) 

where  u  =  2nf,  f  is  the  frequency  of  the  wave,  uq  is  the  permeability  in 
vacuum,  and  a  sinusoidal  time  dependence  is  assumed.  In  the  i-th  layer  of  the 
atmosphere,  n(x,z)  will  satisfy  the  inhomogeneous  Helmholtz  equation: 


2-3 


ECAC-TR-88-001 


Section  2 


2  2  2 
v  n.  +  km.  (z)  n. 

1X1 


J(x,z)/( juy  ) ,  1  <  i  <  L 


(2-5) 


where  the  source  term  is  given  by 


J(x,z)  =  -  j  oj  uQ  p  6 ( x-x '  )  6(z-zt)  (2-6) 

1  /2 

k  =  u)(uqEq)  is  the  wave  number  in  free  space,  and  eq  is  the  permittivity  in 
vacuum. 

The  ni  satisfy  the  following  boundary  conditions  which,  using  Equation 
2-4,  expresses  the  continuity  across  the  region  boundaries  of  the  components 
of  £  and  H  parallel  to  the  boundaries: 

n.  =  n.  . 

i  1+1 

for  z  =  z .  ,  1  <  i  <  L- 1  (2-7 ) 

i 

an./az  =  an.  ,/az 
i  i+l 

Since  it  is  assumed  that  the  ground  conductivity  is  infinite,  the  boundary 
condition  at  the  ground  is 


n  x  lx  =  0 


for  z  =  f(x) 


(2-8) 


where  n  is  the  unit  outward  normal  to  the  surface  and  f(x)  is  a  periodic 
function  that  describes  the  earth's  rough  boundary.  Using  Equation  2-4,  this 
can  be  expressed  as 
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n^(x,z)  =  0  for  z  =  f(x) 


(2-9) 


When  z  is  small,  n^(x,z)  may  be  expanded  in  a  Taylor  series  about  z  =  0  to 
yield 


n.(x,0)  +  [an  (x,o)/az]  f(x)  +  [a2n1(x,0)/az2]  f(x)2/2  +  ...  =  0 

(2-10) 


which,  together  with  Equation  2-7,  are  the  boundary  conditions  for  the 
problem. 

PLANE  WAVE  SPECTRUM  REPRESENTATION 


For  the  smooth  earth  problem,  a  Fourier  transform  formalism  may  be  used 
to  eliminate  the  horizontal  x  variable,  thereby  transforming  the  partial 
differential  equation  into  an  ordinary  differential  equation  (Reference 
1-4).  This  formalism  requires  horizontal  homogeneity  for  its  implementation, 
which  is  not  present  in  the  problem  under  consideration  because  of  the 
presence  of  ground  roughness.  (This  is  not  to  be  confused  with  the 
atmospheric  layers,  which  are  horizontally  homogeneous  and  which  thereby 
permit  the  use  of  the  Hertz  potential  formalism.)  Therefore,  an  analogous 
method  will  be  used  that  is  valid  for  the  periodic  type  of  boundary  that  is 
assumed.  This  method  will  reduce  to  the  Fourier  transform  formalism  in  the 
zero  roughness  limit. 

The  source  function  J(x,z)  may  be  written  as  a  sum  of  its  plane  wave 
spectral  components.  Since  these  components  are  continuously  distributed, 
this  sum  is  expressed  as  an  integral: 
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CD 

J(x,z)  =  — —  f  du  e^uX  J(p,z)  (2-11) 

(2lt)  -0° 

where,  using  the  definition  of  the  Fourier  transform  and  Equation  2-6, 


00 

J  (p,z)  =  — —  f  dx  e”j;jX  J(x,z) 

(2tt)  -® 


=  -  j  uj  pQ  p  6(z 


zT) 


-JPX 


(2-12) 


Each  spectral  component  3(  p ,  z  )e-^x/(2n )  may  be  considered  an  independent 

source  contributing  to  the  electric  field.  Let  the  solution  of  the  problem 

'v  iu(  x-x 1  ) 

for  this  component  of  the  source  be  given  by  lK(p,x,z)eJ  / ( 2tt )  in  the 

i-th  atmosphere  layer;  then,  by  superposition,  the  solution  of  the  problem  for 
all  the  source  components  is 


n.(x,z)  =  — —  J  dp  e^x  x  ^  5.  (u,x,z) 

(2it)  — 


(2-13) 


This  equation  differs  from  the  familiar  Fourier  transform  in  that  IK  is  a 
function  of  x  as  well  as  p  and  z. 

iy(  X-X  *  ) 

Since  the  solution  of  the  problem  in  the  i-th  layer  is  IK(u,x,z)e 

when  the  source  function  is  3(p,z)e-'uX,  these  may  be  substituted  into  Equation 

2-5  for  the  n.  and  J,  respectively.  Carrying  out  the  differentiation  of  the 
lux  ^ 

eJ  factor  and  using  Equation  2-12  then  results  in 


22  22-  ~  222- 
(3  /  3x  +  a  /  3z  )n .  +  2jp  3I1./3X  +  (k  nu  -  p  )n. 

=  -  p  6(  z  -  Zj)  (2-14) 
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i  y  (  X  — X  '  ) 

The  eJ  factor  no  longer  appears  in  Equation  2-14.  The  only  source  of 

variation  in  the  x  direction  is  due  to  the  periodic  rough  earth  boundary 
condition.  Hence,  the  variation  of  ff  with  x  is  periodic  and  may  be  written  in 
the  form: 


nL  (y,x,z) 


„  _  (n),  ,  jicnx 

h  Fi  (lJ-z)  e 


(2-15) 


Substituting  this  form  of  the  solution  into  Equation  2-14  results  in 


I  {d2F.^n)/dz2  +  [k2m.2  -  (y  +  <n)2]  F.(n)}  e^<nx 
n  1  i  i  i  1 

=  -p  6  (z-zT)  (2-16) 


This  equation  could  have  been  obtained  directly  from  Equation  2-5  by 
assuming  a  solution  of  the  form: 


r  F.(n)(u,z) 


J(w  +  <n  )x 


(2-17) 


This  equation  is  a  Floquet  series.  It  is  a  discrete  sum  of  plane  waves 
because  the  boundary  conditions  are  periodic. 

Since  the  m^  are  not  functions  of  x,  the  coefficients  of  exp(j<nx)  in 
Equation  2- 1 6  are  independent.  Therefore,  the  sum  of  N  terms  are  M 
independent  differential  equations.  The  result  is 

d2F^nVdz2  +  [k2rtu(z)2  -  i/n^2]  F^  ^  =  -fi^p  ~  ZT^’  n 

(2-18) 
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where 


M 


(n) 


p  +  icn 


(2-19) 


From  Equation  2-17,  it  is 
each  atmospheric  region  except 
conditions  in  Equation  2-7  can 
IK  replaced  by  F^n^: 


seen  that  the  expression  for  IK  is  the  same  in 
for  the  F^n^  factor.  Hence,  the  boundary 
be  rewritten  in  the  same  form  but  with 


r  (n)  (n) 

F .  -  F  .  - 

1  1+1 


3F. (n)/3z 
1 


3F. 


(n) 


if  1 


/3z 


for  z  =  z..,  1  <  i  <  L-1  (2-20) 


The  expression  of  the  boundary  condition  of  Equation  2-10  is  somewhat 
more  involved  because  of  its  dependence  on  x.  Using  Equation  2-17,  this 
boundary  condition  becomes 


Z  eJ<nx  {F1  (v*,0)  f  3F1(n)(y,0)/3z  f(x) 


+  3  F1(n)(y,0)/3z2  f(x)2/2  f  ...} 


(2-21) 


or, 


ignoring  terms  of  order  f  ( x  P  and  higher, 


I  (e-|,cnx  F <">(„, 0)  . 
n  1  1 

-  J  Y<n)2F1(n)(u,0) 


9F, <">(„, o)/az  i  p  ej<(n,*n)x 

1  m  m 

z  s  ej<(m+n)x}  =  0 

m  m  ‘ 


(2-22) 
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where 


(n)2  .2  2  (n)2 

Y  =  k  m  -  u 


(2-23) 


and  Equation  2-16  was  used.  pm  and  sm  are  the  coefficients  of  the  Fourier 

p 

series  of  f(x)  and  f  (x),  respectively,  which  are  small  for  small  roughness. 
Equation  2-22  will  be  simplified  by  using  the  orthonormality  of  functions  of 
the  form  exp(  j2irnx/d) ,  expressed  as: 

1  j  j2,(m-n)x/d  =  (2-2H) 

d  J  mn 

Multiplying  Equation  2-22  by  e  J  (where  r  is  an  integer)  and 
integrating  the  result  over  the  period  d  yields 


F1 (r) ( y ,0)  +  Z  3F1(n)(y,0)/3z  Pr_p 
-  i  r  y  (n)2  F^ (n)  (u,0)  sr_n  =  0,  all  r 


(2-25) 


Therefore,  to  obtain  a  solution  to  the  rough  earth  problem,  Equation  2-18  must 
be  solved  using  the  boundary  conditions  in  Equations  2-20  and  2-25. 

SOLUTIONS  OF  THE  EQUATIONS 


By  substituting  Equation  2-2,  Equation  2-18  can  be  cast  into  the  form  of 
Stokes'  equation,  the  solution  of  which  is 


F  <"> 

i 


(u,z) 


A.(n)K1(q.(n)) 


B.(n)K2(q.(n)) 


FP  6iP5nO, 


all  n 


(2-26) 
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where 


(2-27) 

K1  and  Kp  are  linear  combinations  of  modified  Hankel  functions  of  order  one- 
2- 1 

third,  and  P  is  the  number  of  the  layer  containing  the  source  current.  The 
Ai(n)  and  B^(n)  are  functions  of  the  parameter  y  and  are  found  by  invoking  the 
boundary  conditions.  It  is  to  be  noted  that,  although  y  entex'S  the 
differential  Equation  2-18  in  the  form  y^n^,  A^n^  and  cannot  be 

written  simply  as  functions  of  y^n^  since  they  must  be  found  by  solving  a 
linear  system  of  equations  that  include  the  boundary  condition  of  Equation 
2-25,  which  contains  all  values  of  n. 

The  particular  solution  Fp  of  the  inhomogeneous  form  of  Equation  2-18 
(i.e.,  n  =  0)  is  given  by  either  of  the  following  (Reference  1-3): 


(n) 


q(<n)<Z> 


2/3 


tana . 
i 


2,  .  y 

.m  (z)  -  — 

k 


F 

P 


Bp  K,(qp<)  K^(  ip>) 


-  Rp  K,(qp>)  *2<V’ 


where 


(2-28) 


Rp  =  P/Wqp ' (z ) 


(2-29) 


ip< 


q  (0)(min{z,z  }), 


*P> 


q  ^0)(max{z,z  } ) 


(2-30) 


‘^Harvard  Computational  Laboratory,  Tables  of  Modified  Hankel  Functions  of 
Order  One-Third  and  of  Their  Derivatives,  Harvard  University  Press, 
Cambridge,  MA,  19^5. 
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W  is  the  constant  VJronskian  of  K|  and  K2  defined  by 

W  =  W{KrK2}  =  K^qp)  K2'(qp)  -  i<2(qp)  (qp)  (2-3D 


The  primes  in  Equations  2-29  and  2-31  indicate  differentiation  with  respect  to 
the  argument.  In  particular, 


dq.^n\z)/dz  (- 


k  ,2/3  . 

7 - 1  I  tan  a. 

tana .  ;  1 

1 1 


(2-32) 


where  Equations  2-2  and  2-27  were  used.  It  should  be  noted  that  the  right 
side  of  Equation  2-32  is  independent  of  z,  y,  and  the  index  n. 

To  solve  for  the  field  in  the  atmospheric  layers,  the  coefficients 
of  the  general  solution  in  Equation  2-26  must  be  known.  These  will  now 
be  determined  by  solving  a  set  of  simultaneous  linear  equations  obtained  by 
substituting  Equation  2-26  into  the  boundary  condition  of  Equations  2-20  and 
2-25.  This  will  be  accomplished  here  under  the  assumption  that  the  source  is 
not  in  the  atmospheric  layer  bordering  the  ground,  i.e.,  P  >  1.  This 
assumption  removes  the  necessity  of  considering  the  inhomogeneous  solution  of 
Equation  2-26  in  the  roughness  boundary  condition  of  Equation  2-25.  It  has  no 
limiting  effect  for  any  finite  height  of  the  source,  since  the  lowest  layer 
can  always  be  formally  divided  into  two  separate  layers  at  a  point  below  the 
source  height.  Substituting  Equation  2-26  into  Equation  2-25  yields: 


y  [A,(r-n)  U  <r-n)  .  B, ( r~n ^  V  (r-n)l 
L  1  n  1  n 

n 


0,  all  r 


(2-33) 
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(n)  _  (n),  . 

qij  -  qi  uj> 


(2-35) 


Finally,  substituting  Equation  2-26  into  Equation  2-20  yields 


(n)  An)  (n)  (n)  (n)  (n)  (n)  R(n) 

av  ,  A .  +  a-.  ~.B.  +  a~.  ,A.  ,  +  a0.  0.  0B .  . 

2i,2i-1  l  2i,2i  l  2 i,2i+1  1+1  2i,2i+2  1+1 


where 


a  (n)  fl  (n)  +  a  (n)  B(n)  +  a(n)  A(n) 

2i+1,2i-1  i  2i+1,2i  i  2i+1,2i+1  i+1 


(n)  (n) 

+  a2i+1 ,2i+2Hi+2 


1  <  i  <  L-1 


B(n) 

B2i+1 


(2-36) 


(2-37) 


a  (n) 

2i,2i-1 


a  (n) 
2i,2i 


a  <"> 
2i,2i+1 

a  (n) 
2i,2i+2 


(2-38) 
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a  (n) 

2i+ 1 ,2i- 1 

=  K,’Kin)) 

a  (n) 

2i+1 ,2i 

a  (n) 

2i+1 ,2i+1 

--  -  '  v 

(2-39) 

a  (n) 

2i+1 ,2i+2 

--  -  K2,<<'u!,i)<’i:i 1  <’i' 

e  (") 
62i 

=  fP<2=2i>  Sn0  -  StP> 

B(n) 

62i  +  1 

=  dFp(z=Z.)/dz  6n0  (6i+1)P  -  6ip)  /  q.1 

• 

(2-40) 

and  when  i  =  L  -  1 

,  the  term  is  absent  in  order  to 

satisfy 

the 

radiation  condition. 

Equations  2-33, 

2-36,  and  2-37  may  be  written  in  matrix 

form  as 

a  n  =  6  (2-41) 


where  a  is  a  square  matrix,  and  n  and  6  are  column  vectors.  To  illustrate 
their  form,  consider  the  case  in  which  L  =  2,  P  =  2,  and  n  takes  on  the 

values  -1,  0,  and  1.  By  grouping  together  Equations  2-33,  2-36,  and  2-37  for 
each  value  of  n,  a  may  be  written 
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where 


e2  =  f2(z=z1) 

6^  =  [dF2(z=z1 )/d2]/q1 ' 


and 


(2-44) 


(n) 

1 


(n)  (n) 

0  ’  1 


(2-45) 


Equation  2-41  is  solved  for  any  desired  elements  of  n  by  using  standard 
methods  for  solving  a  system  of  linear  equations.  Thus,  following  the 
definitions  introduced  for  the  smooth-earth  case  (Reference  1-3),  the 
coefficients  may  be  written  as  ratios  of  determinants: 

«1<n)  =  I lTal<">l I  /  INI.  Bi(n>  =  lit  Bi(")|l  /  IMI 

(2-46) 

where  is  the  matrix  obtained  by  replacing  the  column  of  a  containing 

the  coefficient  A^n^  by  the  vector  6,  and  is  the  matrix  obtained  by 

replacing  the  column  of  a  containing  the  coefficient  by  the  vector  6. 

The  |  | t  |  j  are  functions  of  the  source  height  z-p,  while  ||a||  is  not. 

The  above  values  of  A^n^,  can  be  substituted  into  Equation  2-26  to 

obtain  the  F^n^,  which  are  used,  in  turn,  to  obtain  ff  from  Equation  2-15. 
This  result  is  then  used  in  Equation  2-13  to  obtain  the  desired  Hertz 
potential.  In  the  event  that  the  ground  were  smooth,  all  the  F^n^  would 
reduce  to  0  for  all  n  *  0.  Since  only  n  =  0  would  be  included,  Fp^(y,  z) 
would  be  obtained.  Far  from  the  source,  the  expression  for  the  electric  field 
relative  to  its  free  space  value  is  proportional  to  n  and  is  precisely  the 
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same  as  that  for  the  cy 1 indrically  symmetric  smooth-earth  case  (Reference 
1-3). 


From  Equations  2-13,  2-15,  and  2-26, 


n  (x,z) 


2u 


.  jvcnx 
Z  eJ 


K1(gi(n)) 


(n) 


) 


(2-47) 


The  above  infinite  integral  may  be  transformed  into  a  contour  integral  by 
closing  the  contour  below,  as  shown  in  Figure  2-2.  Since  the  poles  u  =  4m 
of  the  integrand  are  expected  to  lie  in  the  lower  half  of  the  u  plane,  the 
above  integral  can  be  meaningful  for  |  x  |  >  ■»  only  if  x  is  taken  to  be 
negative.  The  integrand  can  then  be  written  as  the  sum  of  residues: 


n  .  (x,z)  =  ( -2tt j/2ir )  Z  Res 

l  mm 


(2-48) 


where 


Res 


Res 


ii.  (y  ,x,z) 


X  E  exp( - ju  | x I ) 
mm  m 1  ' 


(2-49) 


X 


m 


I M  |/34 


-1 

4  =  4. 


m 


(2-50) 


E  :  Z  e-j<nlXl  E! 
m  n 


(2-5D 


E 

mn 


Mq'"’) 


( n  1 


4  =  4. 


m 


(2-52) 
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Equation  2-48,  therefore,  can  be  written 


n. (x,z) 


-j  l  x  exp[-ju  |  x  I  ]  l  e_JlcnlxiE 
J  m  m  K  J  ml  1  " 


mn 


-j  Z  X  E  exp[-ju  I x I 
J  m  m  m  m1  1 


(2-53) 


SUMMARY 


In  Equation  2-42,  the  matrix  of  the  characteristic  determinant  has  been 
divided  into  submatrixes,  with  those  containing  the  a^^  elements  lying 
along  the  diagonal.  By  using  the  definitions  in  Equation  2-34,  it  is  seen 
that  the  roughness,  characterized  by  the  roughness  amplitude  e,  appears  only 
in  the  rows  containing  the  ai^n^  and  a^0^  terms.  Since  the  index  n  labels 
a  Floquet  mode,  the  elements  of  the  off-diagonal  submatrixes  ( i . e . ,  the 
elements  Un^r”ri\  Vn^r_n^)  are  all  of  at  least  first-order  in  e  and  serve  to 
couple  together  the  various  Floquet  mode  contributions. 

This  is  seen  by  observing  that,  if  all  elements  not  contained  in  the 
diagonal  submatrixes  were  zero,  the  Floquet  modes  may  be  "decoupled"  from  each 
other.  Each  submatrix  on  the  diagonal  would  then  represent  an  independent 
system  of  equations  to  determine  the  contribution  of  a  Floquet  mode  of  the 
problem.  Because  of  the  form  of  the  free  vector  6  in  Equation  2-43,  the  free 
vector  for  each  submatrix  would  then  be  zero  except  for  the  0-th  Floquet 
mode.  The  corresponding  field  coefficients  (i.e.,  B^n^),  therefore, 

would  all  be  zero  except  for  the  0-th  Floquet  mode,  thereby  leading  to  the 
smooth-earth  solution 


n . (x,z) 


'J  lii  \»  Em0  e*P(-JlJX|l 


(2-54) 
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This  does  not  imply  that  eigenvalues  (or  zeroes  of  ||aj|)  of  only  the  0- 
th  order  submatrix  are  present.  Just  as  there  are  zeroes  of  the  0-th  order 
submatrix,  there  are  zeroes  of  the  other  submatrixes  representing  other 
Floquet  modes.  However,  in  the  limit  of  zero  roughness,  only  those  modes 
representing  the  eigenvalues  of  the  0-th  order  submatrix  are  excited. 

Referring  to  Equation  2-52,  the  xA^n^and  Tg.^n)  would  be  non-zero  only  for  n 
=  0  and  only  for  roots  p  =  um  of  ||a||  corresponding  to  the  roots  of  the  0-th 
order  submatrix.  Therefore,  the  sum  in  Equation  2-54  is  over  only  those 
eigenvalues  of  the  0-th  order  submatrix. 

The  locations  in  the  complex  p-plane  of  eigenvalues  pm  of  each  Floquet 
contribution  are  illustrated  schematically  in  Figure  2-3  for  the  case  of  zero 
roughness.  Since  the  only  difference  in  each  submatrix  is  the  index  n  of 

the  set  of  eigenvalues  for  each  value  of  n  would  be  translated  from  the 
set  of  modes  for  n  =  0  by  an  amount  <n  along  the  real  axis.  There  is  no  way 
of  associating  particular  eigenvalues  with  specific  values  of  n  when  the 
roughness  is  large.  When  the  roughness  is  non-zero  but  small,  it  is  expected 
that  their  location  in  the  complex  plane  differs  only  slightly  from  the  zero 
roughness  case  so  that  identification  between  a  given  set  of  eigenvalues  and  a 
particular  Floquet  series  element  can  still  be  made. 
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SECTION/  3 

EFFECTIVE  REFLECTION  COEFFICIENT  AND  WAVEGUIDE 
MODE  ATTENUATION 


GENERAL 


In  the  previous  section,  solutions  to  the  electromagnetic  field  equations 
for  horizontally  polarized  waves  propagating  in  a  duct  situated  over  a 
perfectly  conducting,  rough  earth  were  obtained.  These  solutions  depend  on 
determining  eigenvalues  that  are  the  roots  of  an  infinite  determinant. 

In  this  section,  a  waveguide  mode  formulation  is  used  to  derive  an 
effective  ground  reflection  coefficient.  This  reflection  coefficient  is  used  in 
the  fundamental  waveguide  mode  equation  to  determine  the  eigenvalues  of  a 
waveguide  with  a  periodically  rough  boundary  as  a  perturbation  of  the  smooth 
boundary  values.  The  roughness  is  found  to  cause  attenuation  in  trapped  modes, 
and  this  attenuation  can  be  expressed  in  a  simple  manner.  It  reduces  to  a 
closed-form  expression  when  the  duct  is  approximated  by  a  bilinear  refractivity 
profile.  This  approach  is  entirely  equivalent  to  that  given  in  the  last  section 
but  provides  greater  physical  insight  into  the  waveguide  phenomenon. 

REFLECTION  COEFFICIENTS 


A  horizontally  homogeneous  atmosphere  is  characterized  by  a  continuous 
piecewise  linear  refractivity  profile  with  L  sections  (refer  to  Figure  3-1). 
Using  Equations  2-1  and  2-2,  the  modified  refractivity  of  the  i-th  section  of 
the  profile  can  be  written  as 

M^(z)  s  |[mi2(z)  -  1]  x  10^  ,  1  <  i  <  L  (3-D 

where  m^(z)  is  the  modified  index  of  refraction.  For  sinusoidal  roughness, 
the  ground  surface  satisfies 

f(x)  =  e  sin  (<x  +  $)  (3-2) 
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where  <  -  2n/d  and  d  is  the  spatial  period  of  the  roughness. 

Assume  that  the  region  below  z  =  0  (see  Figure  3-1)  is  half-space 

consisting  of  material  with  the  same  refraction  index  rriQ  as  the  atmosphere  at 
z  =  0,  and  that  a  plane  wave  characterized  by  exp( jv/n^x)exp(-jY^n ^z)  is 

incident  on  the  z  =  0  interface  from  below,  where 

,<n)  =  (k02  -  u<n)2)i  =  k0  sin  e(n)  (3-3) 

=  u  +  xn,  kQ  =  m^k,  3  i3  *“)le  Srazing  angle,  which  may  be  complex, 

and  n  is  a  positive  or  negative  integer.  Since  an  e-ia)t  time  dependence 
is  assumed,  the  wave  is  propagating  in  the  negative  x  direction, 
and  Re(i/n^)  >  0,  Re(y^)  >  0  .  The  total  field  in  the  region  z  <  0  then 
will  be  characterized  by  exp(jy^x)  [exp(-jy^z)  +  exp(jy^n^z)]  , 

where  R^  is  the  corresponding  "upward-looking"  reflection  coefficient. 

Now  consider  the  case  of  a  plane  wave  incident  on  the  rough  surface  from 
above.  The  region  above  is  now  assumed  to  be  homogeneous  with  refraction 
index  mQ,  and  the  incident  field  is  characterized  by  exp(  jv/n^x)exp(  jy^z) . 
Since  the  ground  is  periodic  with  wave  number  k,  there  will  be  waves  reflected 
in  discrete  directions  (see  Figure  3-2),  so  that  the  reflected  field  above  the 
ground  will  have  the  form: 


l  (-«nm  +  Rnm)  exp( ju(m)x)exp(-jy(m)z)  (3-4) 


where  6nm  is  the  Kronecker  delta  function,  and  the  sum  is  over  all  positive 
and  negative  m.  The  Rnm  are  elements  of  a  "downward-looking"  roughness 
reflection  matrix  in  which  the  element  (n,m)  represents  the  portion  of  the 
field  which,  because  of  the  roughness,  is  reflected  into  a  wave  characterized 
by  when  the  incident  wave  is  characterized  by  v/n^.  By  expanding  the 
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Figure  3-2.  Discrete  spectrum  of  plane  wave  reflection 
from  sinusoidal  rough  surface. 


field  and  the  boundary  in  powers  of  the  roughness  amplitude  c  (Reference  3-1') » 

-v 

it  can  be  shown  that,  when  the  horizontally  polarized  E  vector  is  normal  to 
the  direction  of  roughness  variation, 


n  /  \  .  /  2  (n)  .  (n-1 )  (n+1 ) .  3  \  /  - 3  c  \ 

Rn  n  £  =  /2e  Y  'Y  +  y  )  +  0(e  )  (3-5) 

Rn  n+ 1 { E )  =  _J  e  y(in)  ej0  +  °  (e2)  (3-6) 

Rn,n-1(e)  =  -J  e  Y(n)  +  0(E2)  (3-7) 


R 

n,n-m 


U) 


0(e2) 


(3-8) 


*3  1 

Fessenden,  D.  E.,  "Space  Wave  Field  Produced  by  a  Horizontal  Electric 
Dipole  above  a  Perfectly  Conducting  Sinusoidal  Ocean  Surface,"  Third 
International  Conference  on  Antennas  and  Propagation  ICAP  83.  1983- 
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The  upward-  and  downward-looking  reflection  coefficients  are  used  in  the 
following  subsections  to  obtain  the  modes  of  the  atmospheric  waveguide. 

WAVEGUIDE  MODE  EQUATION 

The  modes  of  a  smooth-earth  atmospheric  waveguide  can  be  thought  of  as 
discrete  grazing  angles  8  ,  or  equivalently  discrete  wave  numbers  y  , 
characterizing  plane  waves  at  a  given  altitude  within  the  waveguide.  Choosing 
that  altitude  as  z  =  0,  each  such  plane  wave  will  have  the  property  that, 
after  a  complete  cycle  of  reflection  from  the  upper  boundary  and  reflection 
from  the  lower  boundary,  its  amplitude  and  phase  at  each  value  of  x  along 
z  =  0  will  be  identical  to  those  of  the  original  wave.  This  may  be  expressed 
as 


R(0,(U  ) 

m 


R( u  )  - 
m 


(3-9) 


where  R^  is  the  upward-looking  reflection  coefficient  and  R  is  the  downward¬ 
looking  reflection  coefficient.  For  a  horizontally  polarized  wave  propagating 
over  perfectly  conducting  ground, 


R(y)  =  -1,  all  y 


(3-10) 


so  that  the  um  satisfy 


R(0)(,  ) 
m 


-1 


(3-11) 


Although  Equation  3-9  was  derived  on  the  basis  of  a  single  cycle  of 
reflections  from  the  upper  and  lower  waveguide  walls,  the  result  would  be 
identical  if  any  number  of  cycles  were  used. 
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If  a  sinusoidal  roughness  is  now  considered,  the  situation  is  more 
complex,  since  a  wave  incident  on  the  rough  surface  from  above  is  reflected 
into  an  infinite  set  of  discrete  directions  (see  Figure  3-2).  Limiting  the 
discussion  to  those  wave  components  within  the  propagation  medium 
characterized  by  the  wave  number  Equation  3-9  can  be  generalized  to  a 

rough-surface  waveguide  as 

S(u  ,6)  a  R(0)  R(y(0),6)  -1=0  (3-12) 

m 

where  6  is  a  small  quantity  characterizing  the  roughness, 

R(p(0),6)  =  -1  +  AR(U(0),S)  (3-13) 

is  an  effective  reflection  coefficient,  and  AR  is  the  reflection  coefficient 
correction  due  to  roughness  for  waves  characterized  by  \i^K 

Some  thought  indicates  that  a  component  of  the  ground-reflected  wave  that 
is  not  in  the  specular  direction,  because  of  the  assumed  waveguide  nature  of 
the  propagating  medium,  will  be  reflected  back  to  the  ground  from  the  upper 
wall  of  the  guide  (see  Figure  3-3).  Upon  re-reflection  from  the  ground, 
another  discrete  spectrum  of  waves  is  obtained,  including  a  component 
characterized  by  the  original  y  =  wave  number.  The  same  process  can 

occur  after  any  number  of  ground  reflections. 


Referring  to  Equations  3-5  through  3-8,  note  that  after  a  single 
reflection  from  the  ground,  the  reflection  coefficient  correction  for  the 
component  of  the  reflected  field  is  Rq  q.  This  corresponds  to  wave  BHE  in 
Figure  3-3.  After  an  additional  cycle,  there  is  a  reflection  coefficient 
correction  contribution  in  the  direction  due  to  a  wave  that  was  first 

reflected  into  a  direction  and  then  reflected  back  into  a  i/ ^  direction 


(such  as  BGCJ  in  Figure  3-3).  This  contribution  is 

Rn  [R(m)  R  n]  =  Rn  R  R(m). 

0,m  m,0  0,m  m,0 
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Figure  3-3.  Mode  scattering  out  of  and  into  the  '  specular  direction. 
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If  the  reflection  coefficient  correction  contributions  are  to  be  of  order 
e  ,  it  is  seen  from  Equations  3-6  through  3-8  that  m  can  take  on  only  the 
values  of  ±1.  With  this  limitation,  the  contribution  to  AR  after  an 
additional  cycle  would  occur  as  a  result  of  the  original  y^^  wave  being 
scattered  into  a  y^  wave  after  the  first  reflection,  scattered  specularly 
into  a  y^  wave  after  the  second  reflection,  and  scattered  back  into  a  y^ 
wave  after  the  third  reflection.  This  contribution  to  the  reflection 
coefficient  correction  is 


u,m 


*  R 


A  R  A 
0,m  m,0 


-R(m)2) 


0( e3 ) , 


±1 


This  corresponds  to  BGCIDK  in  the  figure.  In  general,  after  N  +  1  cycles,  the 
reflection  coefficient  correction  contribution  is 


m 


±1 . 


Summing  all  these  contributions  for  m  =  -1  and  m  =  1  yields 

45  --  So,o  -  v,  «_1.0  I  (-Rl',,)i  -  Vi  \o  l  (-R<1,)i 

1  1  (3-14) 

where  the  sums  are  over  all  integers  i  >  1.  Equation  3-14  can  be  consideied 
an  expansion  of  AR  for  small  values  of  jR^~^|  and  |R^|.  As  R^-^  and  R^ 
increase  in  magnitude,  AR  contains  more  of  an  interaction  between  the  ground 
and  the  propagation  medium.  In  such  a  case,  use  of  Rq  q  alone  (which,  like  AR 
in  Equation  3-1,  describes  only  the  energy  scattered  out  of  the  specular 
direction)  is  insufficient  to  compute  AR. 


Substituting  Equations  3-5  through  3-7  in  Equation  3-14  yields 


AR 


-2y 


(0) 


(-D 


R(-1)  -  1 
R(-1}  ♦  1 


+  Y 


(D 


R(  1 )  -  1 
R(1)  ♦  1 


(3-15) 
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where  the  identity  £  (-R)1  =  -R/( 1+R)  was  used.  From  Equation  3-15  it  is 

i 

p 

clear  that  A R  is  proportional  to  e  .  For  future  convenience,  6  will  be 

p 

defined  as  £  e  .  The  relationships  in  Equations  3-13  and  3-14  may  be  used  in 
Equations  3-12  to  obtain  the  eigenvalues  w  . 

EIGENVALUE  PERTURBATIONS 


A  separate  numerical  search  for  the  eigenvalues,  or  roots,  of  Equation 
3-12  is  required  for  each  value  of  e  or  d.  This  eigenvalue  search  is 
extremely  time  consuming.  A  method  is  now  described  that,  to  the  same  degree 
of  accuracy,  permits  the  rapid  calculation  of  the  u_(e)  for  any  set  of  (e,d), 
once  the  wm(0)  have  been  found. 

For  a  small  roughness  amplitude  e,  the  eigenvalues  can  be  expanded  as 


vmU)  =  u_(0)  +  A  u 

m  mm 


(3-16) 


where 


Aw  =  du  (0)/d6  6 
m  m 


(3-17) 


is  the  amount  by  which  the  m-th  eigenvalue  is  changed  as  a  result  of  the 
roughness  through  the  process  of  scattering  energy  out  of  and  back  into  the 
specular  direction.  Now, 


dS  =  S  du  +  Std6 
w  6 


(3-18) 


3-9 


ECAC-TR-88-001 


Section  3 


where  the  subscripts  y  and  6  indicate  partial  differentiation,  and  the  partial 
derivatives  may  be  evaluated  at  y  =  y  (0),  6  =  0.  Choosing  du  and  d6 

such  that  y  remains  an  eigenvalue  requires  that  dS  =  0.  Therefore, 

du  (0)/d6  =  -S./S  (3-19) 

m  6  u 

Substituting  this  in  Equation  3-17  and  using  Equation  3-12  results  in 

A u  =  -AR(u  (0) , 6 )  /  R  (0)(u  (0))  (3-20) 

m  m  y  m 

where  Equations  3-10  and  3-11  were  used. 

Numerical  calculations  indicate  that  for  the  um  representing  the  least 
attenuated  modes  of  propagation, 

y  /kn  1  (3-21) 

m  0 

and  that 

|R{1)  (u)|  >>  1,  |R(_1)(ym)|  <<  1  (3-22) 

1  m  1  1  m  1 


when 


,c/k0  >  ^  ^  um/k0^ 


10 


-4 


(3-231 


From  Equation  3-15,  this  leads  to 


.r  ,  (0)  (-1)  (1).  2 

AR  =  j  y  (y  -  y  )  e 


(3-24) 


Using  Equation  3-21  if  may  be  shown  that  for  ic/kp  not  very  large, 
(±1) 


k0JU/k0)  ( +2  -  </k0)J 


(3-25) 
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Then,  using  Equation  3-3, 

AR  =  it (2  -  X/d)*  -  j  (2  +  X/d)h  [sin  e(0)  kQ2  (x/d)*e2] 

(3-26) 

where  X  =  2ir/kQ.  The  first  factor  in  brackets  on  the  right-hand  side  is  of 

order  unity  for  all  X/d  <  1.  The  second  factor  represents  the  small  quantity 

in  which  R  is  being  expanded.  Comparing  it  with  the  square  of  the  Rayleigh 

factor  of  Equation  3-1,  it  is  seen  that  sin^8  in  the  latter  is  replaced  here 

1 

by  sin  8  times  (X/d)M  For  most  cases  of  interest,  the  Rayleigh  factor  will 
be  smaller.  It  should  not  be  surprising,  then,  if  common  situations  arise  in 
which  additional  terms  in  the  expansion  of  Equation  3-16  are  required. 

Using  Equation  3-26  in  Equation  3-20  yields 

Aum  =  -2[ (2  -  x/d)2  -  j (2  +  x/d)*][sin  e(0)/Ru(0)]  [kQ2(x/d)*]  6 

(3-27) 


TRAPPED  MODES 


The  ratio  sin  O^VR  ^  in  Equation  3-27  depends  on  the  eigenvalue  um. 
For  trapped  modes,  which  are  identified  by  Marcus  and  Stuart  (Reference  1-4), 


Im(sin  6 


0 


(3-28) 


This  ratio  can  be  obtained  from  the  phase  integral  approximation  (Reference 
1-5)  according  to  which 


R(0)  _  e-jU  -  ?”) 


where 


=  Me) 


2k 


j"  (m(z)2  -  C2]5dz 
0 


(3-29) 


(3-30) 
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C  =  itIq  cose  ( 3-3 1 ) 

and  z  =  5(e)  is  the  root  of  m(z)^  -  =  0  and  is  the  height  at  which  a 

ray  with  grazing  angle  9  at  z  =  0  will  experience  a  turning  point.  From 
Equation  3-11,  the  eigenvalues  e  =  6m  satisfy 

ip ( e )  =  (2m  +  3/2 ) it ,  m  =  integer  (3-3 2) 

From  Equation  3-3,  um  =  l<0  cose  so  that 

R  ^  =  k  aR^/acose  =  -  k’1  cote  aR^/3sin8  (3-33) 

y  0  0 

which  should  be  evaluated  at  0  =  0  .  Equations  3-29  through  3-33  can  be 
evaluated  numerically  for  any  desired  surface  duct  profile,  or  can  be 
evaluated  in  closed  form  for  the  piecewise  linear,  two-layer  profile 
illustrated  in  Figure  3-1.  For  this  case,  and  for  all  m, 

sin  6  /R  ^  =  -  l  j  tan  a.  (3-34) 

my  u  J  1 

thereby  providing  a  closed  form  expression  for  Aym  in  Equation  3-27: 

Aym  =  J  [(2  +  x/d)*  +  j (2  -  x/d)*]  tan  kQ2  (\/d)56 

(3-35) 

From  Equations  3-28  through  3-33,  it  can  be  shown  that  for  trapped  modes, 

Re( R  (0))  «  0  (3-36) 

u 

Using  Equations  3-28  and  3-36  in  Equation  3-27  yields 

Im( Ay  )  =  2  (2  -  A/d)*  [j  sin  9(0)  /  R  (0)]  [k.2(X/d)*l  6, 

m  y  U 

d  >  A/2,  trapped  modes  (3-37) 
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PROPAGATION  LOSS  RATE 


The  contribution  of  each  waveguide  mode  to  the  total  field  has  a 
dependence  on  the  distance  |x|  given  by  exp[-jum(E)  |x|],  where  |x|  can 
increase  to  large  values  (References  1-3  and  1-4).  If  Aum  has  a  non-zero 
imaginary  part,  the  exponential  decay  factor  of  exp[-j  Aum  |x|]  can  produce  a 
large  effect  on  the  field  that  varies  with  distance.  The  change  Lm  in 
propagation  loss  for  the  m-th  mode  over  that  which  is  experienced  over  smooth 
earth  is 


Lm  =  -20  log1Q{exp[ Im(Aum) |x| ]}  =  -8.68589  Im(Aum)  |x|  (dB)  (3-38) 

It  is  convenient  to  define  the  rate  of  loss  of  the  m-th  mode  due  to  roughness 
as 


Ld  =  L  / | x I  =  -8.68589  Im(Ay  )  (dB/km)  (3-39) 

Km  m  '  m 


where  Ay  is  in  units  of  km-^. 

The  fact  that  Aum  is  proportional  to  e  permits  a  definition  of  roughness 
loss  per  unit  distance  per  square  of  roughness  amplitude: 


LSm  =  ,0"6  W‘2 


(3-40) 


where  e  is  in  meters. 

Using  Equation  3-37  in  Equations  3-39  and  3-40  yields 
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L  *  8.68589  (10'6)  k  2sin  6  fx/d  (2  -  x/d)  ]  1/2/[2jR  (0)], 

sm  U  m  y 


(dB/km/m  ),  d  >  x/2,  trapped  modes 


(3-41) 


and  using  Equation  3-34  in  Equation  3-41  yields 
LSm  *  1-0857  (10’6)  ltan  ailk02  <x/d)*  (2  ~  x/d^  (dB/km/m2), 


d  >  X/2,  trapped  modes,  bilinear  profile 


(3-42) 


where  X  and  d  are  in  the  same  units,  and  the  kQ,  u ,  and  tan  a-,  are  in  units  of 
km-1.  LSm  as  given  in  Equation  3-42  is  always  positive,  thereby  indicating 
the  detrimental  effect  that  roughness  has  on  the  ability  of  the  trapped  modes 
to  propagate  unattenuated. 

It  is  interesting  to  note  that  if  the  only  contribution  to  the  reflection 
coefficient  correction  came  from  energy  scattered  out  of  the  specular 
direction,  then  AR  =  Rq  q.  Comparing  Equations  3-5  and  3-24,  this  is 
tantamount  to  replacing  -•/ 1  ^  by  in  Equation  3-24.  Since  only  the  real 

part  of  AR  enters  the  expression  of  L,,  for  trapped  modes,  Equation  3-41  would 
be  obtained  for  this  case  as  well.  It  should  be  emphasized  that  Lgm  is 
identical  for  the  partial  scattering  and  the  full  scattering  cases  only  when 
the  modes  are  fully  trapped.  Furthermore,  though  the  loss  rate  with  distance 
may  be  the  same  for  these  modes,  the  roughness  effect  on  the  fields  would  be 
different  since  AR  also  affects  the  phase  of  exp(- JAum| x | ) . 
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To  demonstrate  the  effect  of  roughness  on  the  modes  of  propagation,  a 
plot  is  provided  in  Figure  3-4  of  Lgm  as  a  function  of  the  roughness  spatial 
period  d  for  the  two  trapped  modes  obtained  using  the  surface  duct 
refractivity  provided  in  Figure  3-1.  Lgm  in  these  plots  was  calculated  from 
Equation  3-42  using  a  frequency  of  6.814  GHz,  corresponding  to  X  =  4.4  cm. 

For  these  trapped  modes,  Lgm  exhibits  relatively  large  values  in  the 
vicinity  of  d  =  X.  The  maximum  value  occurs  at  the  value  of  k  at  which 
dLsm/dic  =  0,  which  from  Equation  3-42,  is  seen  to  occur  at  d  =  X. 

SUMMARY 


A  method  for  obtaining  the  effect  of  sinusoidal  roughness  with  a  small 
amplitude  on  the  waveguide  eigenvalues  of  horizontally  polarized  radio  waves 
propagating  in  a  horizontally  homogeneous  atmosphere  has  been  described.  The 
method  is  based  on  the  fundamental  waveguide  mode  equation  and  utilizes  an 
eigenvalue  perturbation  scheme.  The  eigenvalue  perturbations,  as  well  as 
roughness  loss  rates,  are  expressed  in  terms  of  reflection  coefficients  that 
consider  energy  scattered  out  of  as  well  as  into  the  specular  direction. 

A  closed-form  expression  is  provided  for  the  roughness  attenuation  rate 
of  trapped  modes  in  an  atmosphere  approximated  by  a  bilinear  refractivity 
profile.  This  attenuation  rate  was  found  to  be  largest  when  the  propagation 
wave  length  and  the  surface  roughness  spatial  period  are  close  in  value. 
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GENERAL 


The  general  approach  to  the  solution  of  the  basic  equauor.s  involving 
electromagnetic  wave  propagation  in  a  surface  duct  over  a  sinusoidal  rough 
earth  was  presented  in  Section  2.  This  consisted  of  a  waveguide  mode  method 
formulated  for  determining  the  field  strength  over  a  slightly  rough  earth. 

The  eigenvalues  of  the  waveguide  formulation  were  found  as  roots  of  an 
infinite  determinant  that  reduces  to  the  smooth-earth  eigenvalues  as  the 
ground  roughness  e  approaches  zero. 

In  this  section  the  matrix  is  written  in  a  form  such  that,  when  the 
roughness  e  approaches  zero,  the  determinant  becomes  a  product  of 
subdeterminants,  each  representing  a  Floquet  mode  of  propagation.  As  in  the 
last  section,  the  roughness  is  assumed  to  be  sinusoidal,  given  by 
Equation  3-2. 

For  the  case  of  L  =  2,  and  including  for  demonstration  purposes  only  the 
central  9x9  elements, 
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where 


aM) 
a  11 

a(-1} 

a12 

a{-1} 

a21 

a(-1} 

a22 

a(-1} 

a23 
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22 
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31 
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A 0) 

do 
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(D 
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11 

a12 

(1) 

JD 
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21 

a22 

a23 

(D 

31 

a(1) 

a32 

a(1) 

a33 

(4-1) 


a2i<n>  *  *i««l1(n,>,  ■  Ki’(qU<n>)’  1  *  '■  2; 


(4-2) 


W"'  =  -K2(q2,(n)),  a33("> 


-K2'(q2/n))q2Vq1',  all  n 
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‘u("> 

-  vC"’>  1’  -  i 

,(n)2l, 

e.("> 

1 

=  ~(c2/8)  K.(q10{n)) 

Y(n)2 

e2j*, 

c.(n) 

1 

=  -( e2/8)  K.(q10(n)) 

y(n)2 

e-2j*. 

(4-3) 

d  (n) 

1 

-  !«  K.'(q10(">)  Q,’ 

eJ$f 

b.(n) 

1 

*  H  K.'(q10(n))  q,' 

q-J\ 

i  =  1,2;  all  n . 

The  K i ( q )  and  ^(q)  are  linear  combinations  of  modified  Hankel  functions  of 
order  one-third  (Reference  2-1). 

q.  5  q  (n^(z  )  =  (k/|tan  a  | )2/3[m  2(z  )  -  y^2/k2] 

(4-4) 

v(n)  =  u  +  icn  (4-5) 


(n)2  ,2  2  (n)2 

y  =  k  m  -  4 


(4-6) 


Zj  is  the  upper  boundary  of  the  j-th  atmospheric  layer,  Zq  =  0,  nig  =  m^O), 
and  primes  denote  differentiation  with  respect  to  the  argument.  The 
superscript  (n)  denotes  the  index  of  a  Floquet  mode. 

It  has  been  shown  that  in  practical  problems,  the  "power  sum"  expression 


ni(  x  ,z ;  e )  =  {l  |Xm  exp[-j4m(e)|x|]  £  e'j*nl*l  Emn|2)^ 

m  n 
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often  provides  more  accurate  predictions  of  the  duct  field  relative  to  free 
space  (Reference  1-3). 

EQUIVALENT  ZERO  MODE  MATRIX 

When  e  approaches  0,  the  determinant  |  [ a.  {  (  reduces  to  a  product  of 
subdeterminants,  each  representing  a  Floquet  mode  of  the  propagation,  and  Emn 
=  0  except  for  n  =  0.  The  discussion  below  concentrates  on  the  effect  of  the 
roughness  on  the  0-th  Floquet  series  element  and  on  the  contribution  of  that 
element  to  the  total  field.  That  is,  when  the  roughness  is  small,  the  um(e) 
may  then  be  used  to  obtain  the  EmQ[  um(e)  ,z  ,Zjj  c  ]  and  xm L e ) ; e ]  so  that  the  n 
=  0  contribution  to  IK(x,z;e)  may  be  computed. 

The  are  the  roots  of  ||a||.  Referring  to  the  infinite  order  matrix 
form,  Equations  4-1  and  4-3  may  be  used  to  write,  to  second  order  in  e, 


D(n)U)  = 

•u<n,(*> 

a  (n) 
a21 

a  (n) 
a31 

a)2<n)U, 

•22<B> 

a  (n) 

32 

a  (n) 
a23 

a  (n) 

33 

.  (n) 

b11 

K 

b12 

D  <"»  = 

a 

a  (n) 
a21 

a  (n) 

22 

a  (n) 

23 

a  (n) 
a31 

a  (n) 

32 

a  (n> 
a33 

(4-9) 


(4-10) 
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bn(n)  =  K1*(q10(n))q1\  b12(n)  =  K2 ' ( q  1 0 ( } }  q  1  '  (4~11} 


and  the  product  and  sum  in  Equation  4-8  are  over  all  n.  Using  the  first 
expression  in  Equation  4-3, 


[n  D  n  ( e ) 

1  n 


n  d  n  (o)  [i 

in  1 


(4-12) 


so  that  Equation  4-8  may  be  written 


INI  =  n  s< 


(4-13) 


where 


=  (1  -  «  r(n)2|  {D(n)  -ID1"1  |D  <«*'>/0<n*1)]}i  n  >  0 


(4-14) 


,(0)  _  r,  ,  (0)2,  fnf0)  .  n  (0)  fn  (1)/n(1)  .  n  (-D/n(-D 


1-6 


{D'u;  -  6  D  D  vu/Dvu  +  D 

1  a  a  a 


(4-15) 
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;(n)  =  Jl  -  6  r(n)2]  {D(n)  -  6  D  (n)[D  <n"1)/D(n"1)]},  n  <  0 

cL  ct 


(4-16) 


and  all  the  determinants  are  evaluated  at  e  =  0.  Equations  4-14  through 

4-16  are  valid  as  long  as  the  denominators  are  finite. 

The  eigenvalues  corresponding  to  the  n  =  0  Floquet  mode,  therefore,  will 
be  the  roots  y  =  ym(e)  of  S^.  But  determination  of  these  roots  is 
equivalent  to  determining  the  roots  of  the  determinant  of  a  matrix  similar  in 
form  to  that  of  Equation  4-9  with  the  substitution 


au  -  K1 '  -  (j  r  K./q’)  (4-17) 

Y  =  { 6  j  [ D  (1)/D(1)  +  D  ('1)/D('1)]}_1  (4-18) 

where  multiplication  by  the  constant[  1  -  6  y(0)2j-1  was  permitted  under  the 
assumption  that  the  source  is  not  located  in  the  lowest  atmospheric  layer. 

The  matrix  in  Equation  4-9  with  the  substitution  in  Equation  4-17  will  be 
referred  to  as  the  "equivalent  matrix"  for  propagation  over  rough  earth  and  is 
identical  in  form  to  the  characteristic  matrix  for  smooth-earth  propagation 
(Reference  1-3).  y  can  then  be  interpreted  as  the  effective  propagation 
constant  in  the  ground,  and  the  effective  reflection  coefficient  is  given  by 


R 


Y  )  /  (  Y 


(0) 


+  y) 


-1  +  AR 


(4-19) 
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where 


AR  =  2  j  { 


.<°>  (D  .  D  ('')/0('1)l 

a  a 


(4-20) 


Since  Equation  4-19  represents  an  expansion  of  R  about  R  =  -1,  it  is  expected 
to  be  valid  when  AR  is  small  relative  to  unity. 

The  equivalent  matrix  can  be  used  to  determine  all  field  contributions  of 
the  0-th  order  Floquet  mode.  It  cannot  be  used  to  obtain  the  contributions 
from  other  Floquet  modes  that  correspond  to  the  0-th  mode  eigenvalues;  this 
will  remain  a  goal  of  future  study. 

Using  a  rough-earth  version  of  the  fundamental  waveguide  mode  equation, 
it  has  been  shown  that 


AR 


U(‘u 


+  y 


(1) 


2 

c 

2 

(4-21) 


where  R^n^  is  the  reflection  coefficient  "looking  upward"  from  below  the  rough 
surface.  That  is,  if  the  region  below  z  =  0  (see  Figure  4-1)  is  replaced  by 
material  with  the  same  refraction  index  m0  as  the  atmosphere  z  =  0,  and  if  a 
plane  wave  characterized  by  exp(jii  x)exp(-jy  )is  incident  on  the  z  =  0 
interface  from  below,  then  the  total  field  in  the  region  z  <  0  will  be  given 
by  exp(ji/n^x)  [expt-jy^z)  +  R^n^exp(  jy^z)  ] .  Equations  4-20  and  4-21  can 
be  shown  to  be  identical  by  noting  that 
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-  j  y 


-  j  Y 


(n) 


,(n) 


,(n) 


(0)  -  Da(n) 
(0)  +  D  (n) 

a 


(4-22) 


which  can  be  obtained  by  solving  a  linear  system  of  equations  expressing  the 
continuity  of  the  Hertz  potential  and  its  normal  derivative  at  the  ground  and 
at  each  atmospheric  layer  interface. 

The  reflection  coefficient  correction  in  Equation  4-21  includes  energy 
scattered  out  of,  as  well  as  into,  the  specular  direction.  If  only  energy 
scattered  out  of  the  specular  direction  were  included,  then 


AR  =  ^  +  y ^  ^ )  e“V2  (partial  scattering)  (4-23) 

EIGENVALUE  PERTURBATIONS  AND  PROPAGATION  LOSS  RATES 

If  the  above  equivalent  matrix  is  used,  a  separate  eigenvalue  search  is 
required  for  each  value  of  e  or  d.  This  eigenvalue  search  is  extremely  time- 
consuming.  An  eigenvalue  perturbation  method,  therefore,  will  be  utilized 
that,  to  the  same  degree  of  accuracy,  will  permit  the  rapid  calculation  of  the 
um(e)  for  any  set  of  (e,d),  once  the  um(0)  has  been  found. 

For  a  small  roughness  amplitude  e,  the  eigenvalues  of  u  =  um(e), 
corresponding  to  the  0-th  order  Floquet  mode,  satisfy  S^°^(u, 6)  =  0  and  may  be 
expanded  as 


U„U)  = 

m  mm 


(4-24) 
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An 


m 


/S<°> 


V 


6 


(4-25) 


Using  Equation  4-15  with 


D(0)(0)  =  0  when  y  =  y  (0) 

m 


(4-26) 


in  Equation  4-25  yields 


Ay_ 


J(u  ,<)  e* 


m 


(4-27) 


where 


j(vk) 


[D a(°}  /  (4  D,/U;)J  (D 


(0) 


(-1)/D(-1) 


a 


(4-28) 


Each  of  the  determinants  in  this  equation  is  evaluated  at  y  =  ym(0),  6  =  e  =  0 
The  ym(e)  resulting  from  the  use  of  Equation  4-27  with  Equation  4-24  can  then 
be  used  in  the  equivalent  matrix  to  obtain  Xm  and  EmQ  in  Equation  2-53.  Thus, 
once  the  smooth-earth  eigenvalues  have  been  found,  the  field  contributions  of 
the  0-th  Floquet  mode  can  be  determined  for  any  combination  of  e  and  d  without 
performing  a  time-consuming  numerical  search  for  the  new  eigenvalues.  The 
accuracy  of  the  eigenvalues  obtained  using  this  perturbation  method  is  the 
same  as  that  obtained  using  an  eigenvalue  seari  t  of  the  determinant  of  the 
equivalent  matrix,  since  this  determinant  is  itself  the  result  of  an  expansion 
in  the  same  small  quantity. 

Using  the  fundamental  waveguide  mode  equation,  it  was  shown  that 
(Reference  1-10) 
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Ay  =  -  AR  /  R(0) 

m  y 


(4-29) 


where  the  subscript  u  indicates  partial  differentiation  with  respect  to  y. 
From  Equations  4-22  and  4-26,  it  may  be  seen  that,  when  y  =  ym(0), 


R(0)  =  -  2  j  y(0)  D(0)(0)  /  D(0) 

y  y  3. 


(4-30) 


The  consistency  of  Equations  4-29  and  4-27  is  demonstrated  by  substituting 
Equations  4-20  and  4-30  into  Equation  4-29.  Thus,  the  eigenvalue  perturbation 
form  of  Equation  4-27  considers  scattering  both  out  of  and  into  the  specular 
direction.  If  only  scattering  out  of  the  specular  direction  were  to  be 
considered,  Equations  4-23  and  4-30  may  be  used  in  Equation  4-29  to  yield 


Ay  =  j  (D^°\o)  /  4D^]  +  Y^1^)  e2  (partial  scattering) 

m  a  y 

(4-31) 

The  rate  of  loss  due  to  roughness  of  the  m-th  mode  along  the  propagation 
direction  is  defined  by 


L„  =  -8.68589  Im( Ay  )  (dB/km) 

Km  m 


(4-32) 


where  Ay  is  in  units  of  km-^.  The  fact  that  Aym  is  proportional  to  t ^  permits 
a  definition  of  roughness  loss  per  unit  distance  per  square  roughness 
amplitude : 


L<>m  is  particularly  significant  for  trapped  modes  in  a  surface  duct.  For 
propagation  over  smooth  earth,  these  modes  experience  little  or  no  attenuation 
in  the  propagation  direction,  which  accounts  for  the  field  enhancement  within 
the  duct  (Reference  1-4).  For  trapped  modes  within  a  surface  duct 
characterized  by  a  bilinear  refractivity  profile, 

LSm  *  1-0857  (10~6)  | tan  c^j  kQ2  (x/d)*  (2  -  x/d)*  (dB/km/m2),  d  >  X/2 

(4-35) 

where  \  and  d  are  in  the  same  units,  and  kg  =  mgk  and  tan  a1  are  in  units 
of  km-^.  Lgm,  as  given  in  Equation  4-35,  is  always  positive.  This 
demonstrates  that  the  roughness  causes  attenuation  of  the  trapped  modes. 
Equation  4-35  is  valid  for  both  the  complete  scattering  (see  Equation  4-27) 
and  the  partial  scattering  (see  Equation  4-3))  situations. 
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SECTION  5 

DISCUSSION  AND  RESULTS 

To  illustrate  the  effects  of  the  roughness  on  the  modes  of  propagation, 
plots  are  provided  in  Figure  5-1  of  Lgm  as  a  function  of  the  roughness  spatial 
period  d  for  the  two  trapped  modes  that  are  obtained  when  the  surface  duct 
refractivity  profile  of  Figure  5-2  is  used.  LSm  in  these  plots  is  calculated  by 
two  methods:  the  eigenvalue  perturbation  method  using  Equation  4-33  and  the 
equivalent  determinant  method  using  Equation  4-34.  A  frequency  of  6.814  GHz  is 
used,  corresponding  to  X  =  4.4  cm.  This  figure  shows  that  the  eigenvalue 
perturbation  results  for  the  two  trapped  modes  are  identical,  consistent  with 
the  usage  of  Equation  4-35,  which  is  independent  of  the  index  m  and  which  was 
originally  derived  from  eigenvalue  perturbation  results.  Indeed,  the  eigenvalue 
perturbation  plot  for  the  trapped  modes  is  virtually  identical  with  the  closed 
form  prediction  of  Equation  4-35. 

It  is  emphasized  that  the  trapped  mode  results  shown  in  Figure  5-1  are 
unique  to  a  surface  duct.  Values  of  LSm  for  trapped  modes  in  an  elevated  duct 
were  found  to  be  several  orders  of  magnitude  lower,  as  would  be  expected  since 
the  ground  roughness  would  not  influence  trapped  mode  propagation  in  an  elevated 
duct. 


A  discrepancy  is  obvious  in  Figure  5-1  between  the  eigenvalue  perturbation 
and  the  equivalent  matrix  results  for  the  trapped  modes  as  d/x  decreases.  This 
corresponds  to  an  increase  in  the  size  of  the  expansion  parameter  AR,  as  noted 
in  the  discussion  following  Equation  4-20.  Although  the  two  methods  provide  the 
same  degree  of  accuracy,  both  are  valid  only  for  small  values  of  the  expansion 
parameter.  Beyond  this  region  of  validity,  there  is  no  reason  to  expect  the 
results  of  each  method  to  be  similar.  The  closeness  of  AR  to  zero,  or  the 
closeness  of  R  to  -1,  thus  provides  an  indication  of  the  accuracy  of  the 
expansion.  Figure  5-3  illustrates  the  complex  value  of  R  for  the  lowest  order 
trapped  mode,  for  e  =  0.1  m  and  various  values  of  d/X.  The  correlation  between 
the  distance  of  the  R  values  from  -1  and  the  divergence  of  the  results  of  the 
two  computational  methods  is  obvious. 
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Figure 


Re  (R) 


5-3.  The  effective  reflection  coefficient  of  sinusoidal 

ground  for  amplitude  e  =  0.1  m,  frequency  =  6.814  GHz 
for  different  values  of  d/i. 
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Predictions  of  the  field  relative  to  free  space,  using  the  four  lowest 
order  modes,  are  shown  in  Figure  5-4  for  a  6.814  GHz  radio  wave  propagating 
within  the  duct  defined  by  the  refractivity  profile  of  Figure  5-2.  They  are 
illustrated  as  a  function  of  observer  height  at  a  distance  66.7  km  from  the 
source  for  propagation  over  smooth  earth  and  for  three  methods  of  computation 
for  propagation  over  the  sinusoidal  rough  earth:  (1)  the  equivalent  matrix 
method  using  the  entire  reflection  coefficient  correction,  (2)  the  eigenvalue 
perturbation  method  using  the  entire  reflection  coefficient  correction  (i.e., 
Equation  4-27),  and  (3)  the  eigenvalue  perturbation  method  using  the  partial 
reflection  coefficient  (i.e.,  using  Equation  4-31).  The  rough  earth  is 
characterized  by  e  =  0.1  m. 

For  the  predictions  in  Figure  5-4a,  the  source  height  z^  =  29  m  and 
d  =  10  m,  while  in  Figure  5-4b,  d  =  1  m  with  the  same  source  height.  In  both 
figures,  there  is  good  agreement  between  the  predictions  obtained  using 
methods  (1)  and  (2).  The  discrepancy  between  these  predictions  and  the 
partial  reflection  predictions  is  greater  for  the  smaller  value  of  d.  Since 
the  source  is  well  above  the  duct,  only  the  leaky  mode  contributes  to  the 
field.  This  mode  is  the  third  of  the  lowest  order  modes,  thereby  explaining 
the  three  peaks  in  the  plots. 

In  Figure  5-4c,  Zj  =  5  m  and  d  =  10  m.  Since  the  entire  contribution  is 
now  due  to  trapped  modes  that  ate  the  two  lowest  order  modes,  only  two  peaks 
are  apparent  in  the  plots.  Here,  there  does  not  appear  to  be  good  agreement 
between  any  of  the  rough-earth  results.  Figure  5-4d  is  the  comparison  of  the 
power  sum  results  (Equation  4-7).  Agreement  in  this  figure,  as  well  as  the 
fact  that  the  Lgm  in  Figure  5-1  for  each  method  is  virtually  identical, 
indicates  that  a  major  source  of  the  problem  is  due  to  the  relative  phase  of 
each  mode.  When  the  modal  contributions  are  added  without  considering  their 
phase,  the  agreement  between  the  two  complete  reflection  coefficient  results 
improves  except  near  the  null. 
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FIELD  RELATIVE  TO  FREE  SPACE  (UB) 

(b> 


Figure  5-4.  Field  relative  to  free  space  within  duct  of  Figure  5-2 
f  =  6.814  GHz,  |  x  j  =  66.7  km,  over  smooth  earth  and 

over  sinusoidal  rough  earth,  e  =  0.1  m.  (Page  1  of  2). 
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SECTION  6 
SUMMARY 

A  method  has  been  described  for  obtaining  the  effect  of  small  sinusoidal 
roughness  on  horizontally  polarized  waves  that  are  propagating  in  a  surface 
duct  situated  over  large  distances  in  a  horizontally  homogeneous  atmosphere. 
The  method  includes  energy  scattered  out  of,  as  well  as  into,  the  specular 
direction  of  modes  of  the  corresponding  smooth-earth  waveguide. 

It  was  shown  that  an  eigenvalue  perturbation  scheme  is  as  accurate  as 
directly  obtaining  the  eigenvalues  from  the  equivalent  matrix  but  requires 
much  less  computation  time.  The  relationship  between  the  two  methods  has  been 
detailed,  and  the  eigenvalue  perturbations,  as  well  as  roughness  loss  rates, 
have  been  shown  to  be  consistent  with  results  obtained  using  a  rough-earth 
form  of  the  fundamental  waveguide  mode  equation. 

The  methods  documented  in  this  report  are  a  first  step  toward  developing 
a  mathematical  model  of  propagation  under  the  given  conditions.  They  can  be 
considered  a  research  model  applicable  to  a  limited  set  of  circumstances. 
Further  work  is  needed  to  extend  the  model  and  make  it  more  widely  applicable. 


6-1/6-2 


ECAC-TR-88-001 


Section  5 


LIST  OF  REFERENCES 


1-1.  Kerr,  D.  E.,  "Propagation  of  Short  Radio  Waves,"  MIT  Radiation 
Laboratory  Series,  Vol.  13,  McGraw-Hill,  New  York,  1951. 

1-2.  Rotheram,  S.,  "Radiowave  Propagation  in  the  Evaporation  Duct,"  The 
Marconi  Review."  Vol.  37,  No.  192,  1974,  pp.  18-40. 

1-3.  Marcus,  S.  W.,  "A  Model  to  Calculate  EM  Fields  in  Tropospheric  Duct 
Environments  at  Frequencies  through  SHF,"  Radio  Science,  Vol.  17,  No. 
5,  1982,  pp.  885-901. 

1-4.  Marcus,  S.  W.,  and  Stuart,  W.  D.,  A  Model  to  Calculate  EM  Fields  in 
Tropospheric  Duct  Environments  at  Frequencies  through  SHF,  ESD 
Technical  Report  81-102,  AD-A107710,  DoD  ECAC,  Annapolis,  MD, 

September  1981. 

1-5.  Budden,  K.  G.,  The  Wave  Guide  Mode  Theory  of  Wave  Propagation,  Logos 
Press,  Ltd.,  London,  1961. 

1-6.  Hitney,  H. ,  et  al.,  "Tropospheric  Radio  Propagation  Assessment," 
Proceedings  of  the  IEEE,  Vol.  73,  No.  2,  1985,  pp.  265-283. 

1-7.  Rice,  S.  0.,  "Reflection  of  Electromagnetic  Waves  from  Slightly  Rough 
Surfaces,"  in  The  Theory  of  Electromagnetic  Waves,  M.  Kline,  ed., 
Interscience  Publishers,  New  York,  NY,  1951,  p.  351. 

1-8.  Ament,  W.  D.,  "Toward  a  Theory  of  Reflection  by  a  Rough  Surface," 
Proceedings  of  the  IRE,  Vol.  41,  No.  1,  1953,  pp.  142-146. 

1-9.  Beard,  C.  I.,  "Coherent  and  Incoherent  Scattering  of  Microwaves  from 
the  Ocean,"  IRE  Transactions  on  Antennas  and  Propagation,  Vol.  AP-9, 
No.  5,  1961,  pp.  470-483- 

1- 10.  Marcus,  S.  W.,  "Propagation  in  a  Surface  Duct  Over  a  Two-Dimensional 

Sinusoidal  Earth,"  Radio  Science,  Vol.  23,  No.  6,  1988,  pp.  1039-1047. 

2- 1  Harvard  Computational  Laboratory,  Tables  of  Modified  Hankel  Functions 

of  Order  One-Third  and  of  Their  Derivatives,  Harvard  University  Press, 
Cambridge,  MA,  1945. 

3- 1.  Fessenden,  D.  E.,  "Space  Wave  Field  Produced  by  a  Horizontal  Electric 

Dipole  above  a  Perfectly  Conducting  Sinusoidal  Ocean  Surface,"  Third 
International  Conference  on  Antennas  and  Propagation  ICAP  83,  1983. 


R-1/R-2 


ECAC-TR-88-001 


DISTRIBUTION  LIST  FOR 

AN  ANALYSIS  OF  PROPAGATION  IN  A  SURFACE  DUCT 
OVER  A  PERIODICALLY  ROUGH  EARTH 
ECAC-TR-88-001 
(Page  1  of  3) 


External  No.  of  Copies 

Superintendent  1 

Naval  Postgraduate  School 
Monterey,  CA  93940 

Commander  1 

Naval  Ocean  Systems  Center 

Attn:  Dr.  J.  Rockway  (Code  8105) 

San  Diego,  CA  92152 

Commander  1 

Naval  Ocean  Systems  Center 
Attn:  H.  V.  Hitney 
San  Diego,  Ca  y2152 

M.  R.  Dresp  1 

MITRE  Corp. 

P.0.  Box  208 
Bedford,  MA  01730 

John  L.  Skillman  1 

NSA/R521 

Fort  Meade,  MD  20755 

Douglas  R.  Woods  1 

NSA/R521 

Fort  Meade,  MD  20755 

Frank  T.  Wu  1 

Michaelson  Laboratory 
Naval  Weapons  Center 
China  Lake,  CA  93555 

Perry  Snyder 
Code  5325 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 


1 


DISTRIBUTION  LIST  FOR 
ECAC-TR-;88-001 
(Page  2  Of  3) 

External 

Prof.  James  R.  Wait 

Department  of  Electrical  Engineering 

The  University  of  Arizona 

Tucson,  AZ  85721 

Prof.  L.  B.  Felsen 

Department  of  Electrical  Engineering 
Polytechnic  Institute  of  New  York 
Route  113,  Farmindale,  NY  11735 

R.  A.  Pappert 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 

Prof.  Ezekial  Bahar 
Electrical  Engineering  Department 
University  of  Nebraska 
Lincoln,  NB  68588 

Jeffrey  Knorr 

Naval  Post  Graduate  School 
Monte. 'ey,  CA  93940 

Dr.  Russell  Clarke  DSTO 
Radicwave  Propagation 
Electronics  Research  Laboratory 
P.0.  ->ox  1600,  Salisbury 
South  Australia  5108 

Harr}  Ng 
FSSM 

Natir  ial  Telecommunications  and  Information  Administration 
179  A  imiral  Cochrane  Dr. 

AnnapDlis,  MD  21401 

W.  Fr  .zier 
FSSM 

National  Telecommunications  and  Information  Administration 
179  Admiral  Cochrane  Dr. 

Annapolis,  MD  21401 


ECAC-TR-88-00 1 


DISTRIBUTION  LIST  FOR 
ECAC-TR-88-001 
(Page  3  of  3) 


External  No.  of  Copies 

Defense  Technical  Information  Center  12 

Cameron  Station 
Alexandria,  VA  22314 

Dr.  L.  B.  Wetzel,  Sr.  Scientist  2 

Code  5303 

Naval  Research  Laboratories 
Washington,  DC  20375-5000 

Internal 


XMTC/T .  Treadway 
D0/S.  Cameron 
DQI/R.  Crawford 
DQT/R.  Whiteman 
DC/K.  Kontson 
DCB/R.  Meidenbauer 
DCC/D.  Lowe 
DCF/G.  Romanowski 
DF/S.  Kennison 
DFA/S.  Sorci 
DFB/H.  Knowlton 
DFI/J.  Borders 
DFS/S.  Timerman 
DJA/J.  Smithmyer 
DN/D.  Baran 
DNA/N.  Williams 
DNM/P .  Sautter 
DNN/G.  King 
DR/A.  Fitch 
DRC/B.  Lindsey 
DRD/H.  Riggins 
DRD/M.  Haberman 
DRD/D.  Eppink 
DRD/S.  Marcus 
DRD/W.  Stuart 
DRI/M.  McGarry 
DRS/D.  Williamson 
DRS/P.  Masters 
DS/C.  Gettier 
DSA/P.  Kubaryk 
DSB/J.  Copeland 
DST/J.  Benischek 
DQL 
DQL 


5 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

10 

Camera-Ready 


